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The null surfaces of a spacetime act as one-way membranes and can block information for a 
corresponding family of observers (time-like curves). Since lack of information can be related to 
entropy, this suggests the possibility of assigning an entropy to the null surfaces of a spacetime. We 
motivate and introduce such an entropy functional for any vector field in terms of a fourth-rank 
divergence free tensor P^t with the symmetries of the curvature tensor. Extremising this entropy 
for all the null surfaces then leads to equations for the background metric of the spacetime. When 
P£b is constructed from the metric alone, these equations are identical to Einstein's equations with 
an undetermined cosmological constant (which arises as an integration constant). More generally, 
if P£b is allowed to depend on both metric and curvature in a polynomial form, one recovers the 
Lanczos-Lovelock gravity. In all these cases: (a) We only need to extremise the entropy associated 
with the null surfaces; the metric is not a dynamical variable in this approach, (b) The extremal 
value of the entropy agrees with standard results, when evaluated on-shell for a solution admitting 
a horizon. The role of full quantum theory of gravity will be to provide the specific form of P°$ 
which should be used in the entropy functional. With such an interpretation, it seems reasonable 
to interpret the Lanczos-Lovelock type terms as quantum corrections to classical gravity. 



I. INTRODUCTION 

The strong mathematical resemblance between the dynamics of spacetime horizons and thermodynamics has led 
several authors [l| to argue that a gravitational theory built upon the Principle of Equivalence must be thought of 
as the macroscopic limit of some underlying microscopic theory. In particular, this paradigm envisages gravity as 
analogous to the theory of elasticity of a deformable solid. The unknown, microscopic degrees of freedom of spacetime 
(which should be analogous to the atoms in the case of solids), will play a role only when spacetime is probed at 
Planck scales (which would be analogous to the lattice spacing of a solid [2j). Candidate models for quantum gravity, 
like e.g., string theory, do suggest the existence of such microscopic degrees of freedom for gravity. The usual picture 
of treating the metric as incorporating the dynamical degrees of freedom of the theory is therefore not fundamental 
and the metric must be thought of as a coarse grained description of the spacetime at macroscopic scales (somewhat 
like the density of a solid which has no meaning at atomic scales). 

In such a picture, we expect the microscopic structure of spacetime to manifest itself only at Planck scales or 
near singularities of the classical theory. However, in a manner which is not fully understood, the horizons — which 
block information from certain classes of observers — link certain aspects of microscopic physics with the bulk 
dynamics, just as thermodynamics can provide a link between statistical mechanics and (zero temperature) dynamics 
of a solid. (The reason is probably related to the fact that horizons lead to infinite redshift, which probes virtual high 
energy processes; it is, however, difficult to establish this claim in mathematical terms). It has been known for several 
decades that one can define the thermodynamic quantities entropy S and temperature T for a spacetime horizon [J]. 
If the above picture is correct, then one should be able to link the equations describing bulk spacetime dynamics with 
horizon thermodynamics in a well defined manner. 
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There have been several recent approaches which have attempted to quantify this idea with different levels of success 
P) S @ ■ An explicit example was Q the case of spherically symmetric horizons in four dimensions. In this case, 
Einstein's equations can be interpreted as a thermodynamic relation TdS = dE + PdV arising out of virtual radial 
displacements of the horizon. More recently, it has been shown 0] that this interpretation is not restricted to Einstein's 
General Relativity (GR) alone, but is in fact true for the case of the generalised, higher derivative Lanczos-Lovelock 
gravitational theory in D dimensions as well. Explicit demonstration of this result has also been given for the case of 
Friedmann models in the Lanczos-Lovelock theory \£Q as well as for rotating and time dependent horizons in Einstein's 
theory [loj . In a related development, there have been attempts to interpret other gravitational Lagrangians (like 
f(R) models) in terms of non-equilibrium thermodynamics 

In standard thermodynamics, extremisation of the functional form of the entropy (treated as a function of the 
relevant dynamical variables) leads to the equations governing the equilibrium state of the system. This suggests that 
in the context of gravity as well, one should be able to define an entropy functional, which — on extremisation — 
will lead to the equations describing the macroscopic, long-wavelength behaviour of the system (which in this case is 
the spacetime.) That is, if our analogy is to be taken seriously, we should be able to define an "entropy functional for 
spacetime" , the extremisation of which should lead to the gravitational field equations for the metric of the spacetime. 
At the lowest order, this should give Einstein's equations but the formalism will continue to be valid even in the 
quantum regime. It is important to recall that — even in the case of ordinary matter — there is no such thing as 
'quantum thermodynamics'; only quantum statistical mechanics. What quantum theory does is to modify the form of 
the entropy functional (or some other convenient thermodynamic potential, like free energy); the extremisation now 
leads to equations of motion which incorporate quantum corrections. Similarly, we expect our entropy functional to 
pick up corrections to the lowest order term, thereby leading to corrections to Einstein's equations. In that sense, 
this approach is very general and this is what we will develop in this paper. 

A crucial point to note is that the concept of entropy, both in the standard thermodynamics as well as in the 
gravitational context, stems from the fact that certain degrees of freedom are not observable for certain classes of 
observers. (Throughout the paper, we will use the term observers to mean family of time like curves, without any 
extra connotations.) In the context of a metric theory of gravity this is inevitably linked with the existence of one-way 
membranes which are provided by null surfaces in the spacetime. The classical black hole event horizons, like e.g., the 
one in the Schwarzschild spacetime, are special cases of such one-way membranes. This also leads to the conclusion — 
suggested by several authors e.g., [a, [l2j — that the concept of entropy of spacetime horizons is intrinsically observer 
dependent, since a null surface (one-way membrane) may act as a horizon for a certain class of observers but not 
for some other class of observers. In flat Minkowski spacetime, the light cone at some event can act as a horizon for 
the appropriate class of uniformly accelerated Rindler observers, but not so for inertial observers. Similarly, even the 
black hole horizon (which can be given a 'geometrical' definition, say, in terms of a Penrose diagram) will be viewed 
differently by an observer falling into the black hole compared to another who is orbiting at a radius r > 2M. The 
fact that the dynamics of the spacetime should nevertheless be described in an observer independent manner leads to 
very interesting consequences, as we shall see. 

The above discussion also points to the possibility that the null surfaces of the spacetime should play a key role 
in the extremum principle we want to develop. It is also important, within this conceptual framework, that the 
metric is not a dynamical variable but an emergent, long wavelength concept [l3T |. In this paper, we will construct 
such an entropy functional in a metric theory of gravity and derive the equations resulting from its extremisation. 
We will show that not only Einstein's GR, but also the higher derivative Lanczos-Lovelock theory can be naturally 
incorporated in this framework. The formalism also allows us to write down higher order quantum corrections to 
Einstein's theory in a systematic, algorithmic procedure. 

The paper is organized as follows: In section 2 we motivate a definition for the entropy functional S related to a 
vector field £ a in the spacetime in the context of the Einstein and Lanczos-Lovelock theories, and elaborate upon the 
variational principle we employ to determine the spacetime dynamics from this functional. In section 3 we compute the 
extremised value of S[£] and show that under appropriate circumstances it is identical to the expression for the horizon 
entropy as derived by other authors in the context of Lanczos-Lovelock theories, thus justifying (at least partially) the 
name 'entropy functional'. Section 4 rephrases our results in the language of forms, to give its geometrical meaning. 
We conclude in section 5 by discussing some implications of our results. 

II. AN ENTROPY FUNCTIONAL FOR GRAVITY 

Our key task is to define a suitable entropy functional for the spacetime. Since this is similar to introducing 
the action functional for the theory, it is obvious that we will not be able to derive its form without knowing the 
microscopic theory. So we shall do the next best thing of motivating its choice. (If the reader is unhappy with the 
motivating arguments, (s)he may take the final form of the entropy functional in Eq.([T]) below as the basic postulate 
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of our approach!) 

The first clue comes from the theory of elasticity. We know that, in the theory of elasticity [lj], the key quantity 
is a displacement vector field £ a (x) which describes the elastic displacement of the solid through the equation x a — > 
x a + £ i a (x). (Of course, in elasticity, one usually deals with three-vectors while we need to work in D-dimensions!). All 
thermodynamic potentials, including the entropy of a deformable solid can be written as an integral over a quadratic 
functional of the displacement vector field, which can capture the relevant dynamics in the long-wavelength limit. 
In the context of gravity, the "solid" in question is spacetime itself [l5| . The crucial difference from the theory of 
elasticity is the following: In elasticity, extremising the entropy function will lead to an equation for the displacement 
field and determine £ a . In the case of spacetime, the equations should determine the background metric. This is a 
nontrivial constraint on the structure of the theory and we will show how this can be achieved. 

In the case of elasticity as well as gravity, we will expect the entropy functional to be an integral over a local 
entropy density, so that extensivity on the volume is ensured. In the case of an elastic solid, we expect the entropy 
density to be translationally invariant and hence depend only upon the derivatives of £ a quadratically to the lowest 
order. We would expect this to be true for pure gravity as well, and hence the entropy density should have a form 
Pab Cd ^ c£ a V d£, b ', where the fourth rank (tensorial) object P ab cd is built out of metric and other geometrical quantities 
like curvature tensor of the background spacetime. But in the presence of non-gravitational matter distribution in 
spacetime (which, alas, has no geometric interpretation), one cannot demand translational invariance. Hence, the 
entropy density can have quadratic terms in both the derivatives V a £ fc as well as £ a itself. We will denote the latter 
contribution as T ab ^ a ^ b where the second rank tensor T ab (which is taken to be symmetric, since only the symmetric 
part is relevant to this expression) is determined by matter distribution and will vanish in the absence of matter. 
(We will later see that T ab is just the energy momentum tensor of matter; the notation anticipates this but does not 
demand it at this stage.) So our entropy functional can be written as: 

S[(\ = [ d D x^—g{AP ab cd V c ^V d e~T ab Ce) , (1) 
Jv 

where V is a £>-dimensional region in the spacetime with boundary dV, and we have introduced some additional 
factors and signs in the expression for later convenience. We will now impose two additional conditions on P ab cd and 
T ab . (i) For the case of the elastic solid, the coefficients of the quadratic terms are constants (related to to the bulk 
modulus, the modulus of rigidity and so on). We take the analogues of these constant coefficients to be quantities 
with vanishing covariant divergences. That is, we postulate the "constancy" conditions: 

V b P/ cd = = V a T ab . (2) 

(ii) The second requirement we impose is that the tensor P a bcd should have the algebraic symmetries similar to the 
Ricmann tensor R abc d of the D-dimensional spacetime; viz., P abc d is antisymmetric in ab and cd and symmetric 
under pair exchange. Eq.@ then implies that the P abcd will be divergence- free in all its indices. Because of these 
symmetries, the notation P^ with two upper and two lower indices is unambiguous. 

In summary, we associate with every vector field £ a in the spacetime an entropy functional in Eq.|T]), with the 
conditions: (i) The tensor P a bcd is built from background geometrical variables, like the metric, curvature tensor, etc. 
and has the algebraic index symmetries of the curvature tensor. It is also divergence free, (ii) The tensor T ab is related 
to the matter variables and vanishes in the absence of matter. It also has zero divergence. One key feature of the 
functional in Eq. (fTJ) is that the entropy associated with null vector fields are invariant under the shift T ab — > T ab + pg ab 
where p is a scalar. This fact will play an interesting role later on. 



1. Explicit form of P abcd 

Obviously, the structure of the gravitational sector is encoded in the form of P abcd and we need to consider the 
possible choices for p abcd which determine the form of the entropy functional. In a complete theory, the form of 
pabcd w jjj determined by the long wavelength limit of the microscopic theory just as the elastic constants can 
- in principle — be determined from the microscopic theory of the lattice. However, our situation in gravity is 
similar to that of the physicists of the eighteenth century with respect to solids and — just like them — we need to 
determine the "elastic constants" of spacetime by general considerations. Taking a cue from the standard approaches 
in renormalization group, we expect p abcd to have a derivative expansion in powers of number of derivatives of the 
metric: 

(1) (2) 

P abcd ( 9ij , R ljkl )= Cl P abcd ( 9ij ) + c 2 P abcd ( 9ij ,R ijkl ) + ..., (3) 



4 



where c\,C2,- ■ ■ are coupling constants. The lowest order term must clearly depend only on the metric with no 
derivatives. The next term depends on the metric and curvature tensor. Note that since P abcd is a tensor, its 
expansion in derivatives of the metric necessarily involves the curvature tensor as a "package" comprising of products 
of first derivatives of the metric (the TT terms) and terms linear in the second derivatives (<9F), where F symbolically 
denotes the Christoffcl connection. Higher order terms can involve both higher powers of the curvature tensor, as well 
as its covariant derivatives. 

These terms can, in fact, be listed from the required symmetries of P a bcd- For example, let us consider the possible 
fourth rank tensors p ahcd which (i) have the symmetries of curvature tensor; (ii) are divergence-free; (iii) are made 
from g ab and R a bcd but not derivatives of R a bcd - If we do not use the curvature tensor, then we have just one choice 
made from the metric: 

P a c b d = ^(S a A-S a d S b c ). (4) 

OZ7T 

We have fixed an arbitrary constant in the above expression for later convenience. Next, if we allow for p abcd to 
depend linearly on curvature, then we have the following additional choice of tensor with the required symmetries: 

Pfd = ^ ( R cd - G a J b d + G b c S a d + R a d S b - R b d 6 a c ) . (5) 

We have again chosen an arbitrary constant for convenience, but in this case the constant can always be specified in 
the factor c 2 of Eq.(g]). 

The expressions in Eq.Q and Eq.([5]) can be expressed in a more illuminating form. Note that, the expression in 
Eq.flJ) is just 

(!) 1 1 

where we have the introduced the alternating or 'determinant' tensor S b ^ b \ Similarly, the expression in Eq.Q above 
can be rewritten in the following form: 



( 2 ) 11 

pa\a 2 fttif2°3°4 f>b 3 b 4 fj\ 

r b 1 b 2 -[.67r 2 1 b2 bs bi a3Q4 ' ' 



where we have again introduced the alternating tensor S^f^ 

ca 1 a 2 a 3 a 4 _ ,ci ■■■cp_4«l 020304 , , , , , (R\ 

u b 1 b 2 b 3 b 4 ~ (£)_4)l t ci---c D _ 4 fcib 2 0304 ) \°) 

The alternating tensors are totally antisymmetric in both sets of indices and take values +1, —1 and 0. They can 
be written in any dimension as an appropriate contraction of the Levi-Civita tensor density with itself [l6| . (In 4 
dimensions the expression in Eq.© is essentially the double-dual of R a bcd-) We see a clear pattern emerging from 
Eq.([S]) and Eq.© with the m-th order contribution being a term involving (to — 1) factors of the curvature tensor. 
Following this pattern it is easy to construct the m-th order term which satisfies our constraints. This is unique and 
is given by 



(m) 

P cd rv- X c da,3...a 2m T>b 3 b 4 T>b 2m - 1 b 2m (q\ 

r ab 01 u abb 3 ...b 2m n a 3 a 3 ' n a 2m _ia 2 ™ ' \° J 



These terms have a close relationship with the Lagrangian for Lanczos-Lovelock theory, which is a generalised higher 
derivative theory of gravity. Before proceeding further, we shall briefly recall the properties of Lanczos-Lovelock 
theory and describe this connection. 

The Lanczos-Lovelock Lagrangian is a specific example from a general class of Lagrangians which describe a (possibly 
semiclassical) theory of gravity and are given by 

C = Q bcd R\ cd , (10) 

where Q bcd is the most general fourth rank tensor sharing the algebraic symmetries of the Riemann tensor R a bcd and 
further satisfying the criterion VbQa Cd = (Several general properties of this class of Lagrangians are discussed in 
Ref. [171]). The Z?-dimensional Lanczos-Lovelock Lagrangian is given by [18j a polynomial in the curvature tensor: 

K 1 

r( D ) — \ , n r( D ) ■ r( D ) — 9~ m R a ^ a ^ ■■■"2m yb 1 b 2 r>b 2m ^ 1 b 2m 

L - — CmLj rn i >~m ~ °b 1 b 2 ...b 2m rL a 1 a 2 rL a 2m _ 1 a 2m ! K 13 -) 

m—1 
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where the c m are arbitrary constants and Cm is the m-th order Lanczos-Lovelock term. Here the generalised 
alternating tensor 5"\ is the natural extension of the one defined in Eq.© for 2m indices, and we assume D > 2K + 1. 
The m-th order Lanczos-Lovelock term Cm given in Eq. tfTTj) is a homogeneous function of the Ricmann tensor of 
degree m. For each such term, the tensor Q a bcd defined in Eg. (fit))) carries a label m and becomes 

(m)^") cd _ 1 r,-mxcda 3 ...a 2m pb 3 fc 4 Tfb 2m -lb2m { ~\ 9\ 

b ~ 167T abb 3 — b 2m -"-oaas ^fem-iaam • v ' 

The full tensor Q ab cd is a linear combination of the ^ m ^Q ab cd with the coefficients c m . Einstein's GR is a special case of 
Lanczos-Lovelock gravity in which only the coefficient C\ is non-zero. Since the tensors ^ n 'Q„i ) cd appear linearly in the 
Lanczos-Lovelock Lagrangian and consequently in all other tensors constructed from it, it is sufficient to concentrate 
on the case where a single coefficient c m is non-zero. All the results that follow can be easily extended to the case 
where more than one of the c m are non-zero, by taking suitable linear combinations of the tensors involved. Hence, 
to avoid displaying cumbersome notation and summations, we will now restrict our attention to a single m-th order 
Lanczos-Lovelock term Cm , and will also drop the superscript (m) on the various quantities. Comparing with our 
expression in Eq.® it is clear that P * J can be taken to be proportional to Q a ^ h which can be conveniently expressed 
as a derivative of the Lanczos-Lovelock Lagrangian with respect to the curvature tensor. To be concrete, we shall 
take the m-th order term in Eq.© to be: 

dC (D) 

P a ljk = mQS Jk = M^ k = . (13) 

01t ijk 

This equation defines the divergence-free tensor M a z: > k , where the partial derivatives are taken treating g ab , T a bc and 
cd as independent quantities. The numerical coefficients are chosen for later convenience and can be, of course, 
absorbed into the definitions of the c m . With this choice, we have completely defined the geometrical structure of the 
entropy functional, except for the coupling constants c m which appear at each order [19J. 

Just to see explicitly and in gory detail what we have, let us write down the entropy functional in the absence of 
matter {T ab = 0), correct up to first order in the curvature tensor in p abcd . To this order, our entropy functional ([TJ) 
takes the form S = Si + S2 where 

si® = / v ^((v c o 2 -v a £ fc v b r) 

SaK] = Ca / d D x(RZiV c Z a V d e ~ (G c a + R c a )(V c Z a Vbe - Vct b V b n) (14) 
Jv 

where C2 is a coupling constant and we have used equations Q and ([5|). We will later see that the entropy given by 
Si leads to Einstein's equations in general relativity while S2 and higher order terms can be interpreted as corrections 
to this. If we choose £ a to be the normal vector field of a sequence of hypersurfaces foliating the spacetime, then the 
integrand in Si has the familiar structure, (TrK) 2 — Tr(K 2 ) where K a b is the extrinsic curvature. (This could offer 
an alternative interpretation of ADM action; however, we will not discuss this aspect here except to cast our results 
in the familiar language when appropriate.) The expression in Eq. (|14p can be further simplified by integrating it by 
parts where appropriate and writing the right hand side as a sum of a contribution from the bulk and a surface term. 
The general expression after such a splitting is given later in Eq. (f2U)) and can also be found in Section 4 where it 
arises transparently in the language of forms (see Eq. (|53ll ). 



2. Field equations from extremising the entropy 

Having made these general observations regarding the choice of P"] 1 let us now return to the entropy functional in 
Eq.(p}. This expression is well defined for any displacement vector field £ a . We can, therefore, associate an entropy 
functional with any hypersurface in the spacetime, by choosing the normal to the hypersurface as £ a . Among all 
such hypersurfaces, the null hypersurfaces will play a key role since they act as one-way membranes which block 
information for a specific class of observers. Given this motivation, we will now extremise this S with respect to 
variations of the null vector field £ a and demand that the resulting condition holds for all null vector fields. That 
is, the "equilibrium" configurations of the "spacetime solid" are the ones in which the entropy associated with every 
null vector is extremised. Varying the null vector field £° after adding a Lagrange multiplier A for imposing the null 
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condition £ a (5£ a = 0, we find: 

6S = 2 / d D x^g(4P ab cd V c e (V d 5£ b ) -T ab e5Z b - XgabfS?) 
Jv 

= 2 / d D xV=g~ (4P ab cd V c r (V d <5£ b ) - Taif^) (15) 
Jv 

where we have used the symmetries of P ab cd and T ab and set T ab — T ab + Xg a b- (As we said before such a shift leaves 
entropy associated with null vectors unchanged so the Lagrange multiplier will turn out to be irrelevant; nevertheless, 
we will use T ab for the moment.) An integration by parts and the condition V dP ab cd = 0, leads to 

5S = 2 / d D x^—g (-4P ab cd (V d VcD - T ab £ a ) S^ b + 8 / d^xVh (n d P ab cd (V C C)) 5£ b , (16) 

Jv JdV 

where n a is the D-vector field normal to the boundary dV and h is the determinant of the intrinsic metric on dV. 
As usual, in order for the variational principle to be well defined, we require that the variation S£ a of the null vector 
field should vanish on the boundary. The second term in Eq. lfTBf therefore vanishes, and the condition that S[£] be 
an extremum for arbitrary variations of £ a then becomes 

2P ab cd (V c V d - V rf V c ) e - T ab C = , (17) 

where we used the antisymmetry of P ab cd in its upper two indices to write the first term. The definition of the 
Riemann tensor in terms of the commutator of covariant derivatives reduces the above expression to 

(2P b ^ k R\ 3k -ft)i a = 0, (18) 

and we see that the equations of motion do not contain derivatives with respect to £. This peculiar feature arose 
because of the symmetry requirements we imposed on the tensor P ab cd . We further require that the condition in 
Eq. (|18[) hold for arbitrary null vector fields £ a . A simple argument based on local Lorentz invariance then implies 
that 

2P b ijk R a ijk -T a b =F(gK, (19) 

where F(g) is some scalar functional of the metric and we have absorbed the XS b in T b =T b + X5 b into the definition 
of F. The validity of the result in Eq. (fTT)|) is obvious if we take a dot product of Eq. (fT8")) with £ & . (A formal proof can 
be found in the Appendix (|A ![) .) The scalar F(g) is arbitrary so far and we will now show how it can be determined 
in the physically interesting cases. 

2.1 Lowest order theory: Einstein's equations 

To do this, let us substitute the derivative expansion for p abcd in Eq.Q into Eq. (fl9| . To the lowest order we find 
that the equation reduces to: 

±-Rl - T b a = F{g)8 a b (20) 
sn 

where F is an arbitrary function of the metric. Writing this equation as (G b —%nT b ) — Q(g)5 b with Q = 8ttF— (1/2) R 
and using V a G b = 0, V a T^ = we get d b Q = d b [8irF—(l/2)R] = 0; so that Q is an undetermined integration constant, 
say A, and F must have the form 8irF — (l/2)i? + A. The resulting equation is 

i?£-(l/2)i?,5 b Q = 8^T 6 a + A(5£ (21) 

which leads to Einstein's theory if we identify T ab as the matter energy momentum tensor with a cosmological constant 
appearing as an integration constant. (For the importance of the latter with respect to the cosmological constant 
problem, see Ref. [20j; we will not discuss this issue here.) 

The same procedure works with the first order term in Eq.Q as well and we reproduce the Gauss-Bonnet gravity 
with a cosmological constant. In this sense, we can interpret the first term in the entropy functional in Eq. fH]) 
as the entropy in Einstein's general relativity and the term proportional to C2 as a Gauss-Bonnet correction term. 
Instead of carrying out this analysis explicitly order by order, we shall now describe the most general structure in the 
family of theories starting with Einstein's GR, Gauss-Bonnet gravity etc. — and will show that we reproduce the 
Lanczos-Lovelock theory by our approach. 
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2.2 Higher order corrections: Lanczos-Lovelock gravity 

To see this result, let us briefly recall some aspects of Lanczos-Lovelock theory. It can be shown that (see e.g., [l7| ) 
the equations of motion for a general theory of gravity derived from the Lagrangian in Eq. (|10[) using the standard 
variational principle with g ah as the dynamical variables, are given by 

E ab = ~T ab ; E ab = J^ d 9 gab (V=9£) - 2V m V"Af am » b . (22) 

Here T ab is the energy-momentum tensor for the matter fields. The tensor M abc d defined through M a bcd = (dC / dR a bcd ) 
is a generalisation of the one defined for the Lanczos-Lovelock case in Eq. (fT3"|) . The partial derivatives are as before 

taken treating g ab , T a bc and R a bcd as independent quantities. For the m-th order Lanczos-Lovelock Lagrangian Cm\ 
since M abcd is divergence- free, the expression for the tensor E ab in Eg. ([22]) becomes 

dC (D) 
ab = ~dg^ 2 

where we have used the relation d(^f—g)/dg ab = — (l/2)y/—gg ab . The first term in the expression for E ab in Eq. 
can be simplified to give 



U { ° ] 9a b , (23) 



dC 



(D) 



"'•/.W = M a ^ k R bljk , (24) 



a b 



where the expressions in Eq. ([24j) can be verified by direct computation, or by noting that Cm is a homogeneous 
function of the Riemann tensor R a bcd of degree m. To summarize, the Lanczos-Lovelock field equations are given by 

E a b = \T a b ; E ab = mQ^ k R bijk -^C^g ab . (25) 

Further, diffcomorphism invariance implies that the tensor E ab defined in Eq. ([22|) is divergence-free, ^7 a E b — 0. The 
equations of motion for the matter imply that the energy-momentum tensor T ab is also divergence-free (as required 
by Eq.®). Using these conditions in Eq. (Til)|) together with the choice in Eq. (fr3"|) for P a ljk leads to 

d a F = d a C^ , (26) 

which fixes F(g) as F = Cm + A/87T where A is a constant with the normalisation chosen so as to conform with the 
usual definition of the cosmological constant. The resulting field equations for Lanczos-Lovelock gravity will be: 



16tt 



ijk n« _ }_Xa r (D) 
) ±L ijk 2 6 m 



8nTl (27) 



where we have included a possible cosmological constant, that arises as an undetermined integration constant in the 
defintion of T b a . Taking the trace of this equation, we find that that Cm^ = (2m — D)~ 1 T. In other words, the on-shell 
value of the Lagrangian is proportional to the trace of the stress tensor in all Lanczos-Lovelock theories, just like in 
GR. In the absence of source term, this implies that Cffl — and the equations of motion reduces to P b t ^ k R a i ^ k = 0. 

The case m — 1 with C^J.^ — (l/167r)i? is easily seen to reduce to that of Einstein's gravity. . 

To summarise, if we take the derivative expansion in Eq.([3]) to correspond to a polynomial form in the curvature 
tensor, then it has the form given by Eq. (TT3"|) . In this case, extremising the entropy leads to the Lanczos-Lovelock 
theory. We stress that the resulting field equations have the form of Einstein's equations with higher order corrections. 
In our picture, we consider this as emerging from the form of the entropy functional which has an expansion in powers 
of the curvature. 

Before concluding this section, we want to comment on an interesting property of the entropy functional. The 
derivation of the equations l|19p was based upon a variational principle which closely resembles the usual variational 
principle used in other areas of physics in which some quantity is varied within an integral arbitrarily, except for it 
being fixed at the boundary. Instead of such an arbitrary variation, let us consider a subset of all possible variations 
of the null vector field £ a , given by £, a {x) — ► (1 + e(x))^ a (x)\ namely infinitesimal rescalings of £ a . We assume that 
the scalar e(x) is infinitesimal and also that it vanishes on the boundary dV. In this case it is easy to see that the 
variation of S[£] in Eq. (fT6"|) becomes 

^Ircscalc = 2 J d° R% jk - T?) U^x) (28) 
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Clearly, requiring that the functional S[£] be invariant under rescaling transformations of £ a leads to the same 
requirement as before, namely that Eq. p^)) be satisfied. We can understand the physical motivation behind imposing 
such a symmetry condition on S[S,] as follows. Let us begin by noting the fact that the causal structure of a spacetime, 
which can be thought of as the totality of all possible families of null hypersurfaces in the spacetime, is left invariant 
under rcscalings of the generators of these null hypersurfaces. To see this symmetry, note that any curve in the 
spacetime, with tangent vector field t a = dx a /d\, say, is invariant under the rescaling t a — > f(X)t a of the tangent field, 
where / is some scalar. This is so because a rescaling of the tangent vector field is equivalent to a reparametrization 
of the curve. A null hypersurface TL can be thought of as being 'filled' by null geodesies contained in it, and is hence 
invariant under rescalings of its generator field £°. The result for the full causal structure then follows. The functional 
S[£\ depends only on the generator £ a of some null hypersurface (apart from the metric and matter fields which we 
consider as given quantities). It is therefore natural to demand that this symmetry of the causal structure also be a 
symmetry of S[£\. We do not, however, use this feature in this paper. (For a completely different, purely classical, 
approach to general relativity based on null surfaces, see (HI)- 

III. EVALUATING THE ENTROPY FUNCTIONAL ON-SHELL 

The results in the previous section show that our approach provides - at the least - an alternative variational 
principle to obtain not only Einstein's theory but also Lanczos-Lovelock theory. In the conventional approaches to 
these theories, one can obtain the field equations by varying the action functional. But once the field equations 
are obtained, the extremum value of the action functional is not of much concern. (The only exceptions are in the 
semiclassical limit in which it appears as the phase of the wave function or in some specific Euclidean extension of 
the solution.) In our approach, it is worthwhile to proceed further and ask what this extremal ['on-shell'] value means 
in specific contexts. We will be able to provide an interpretation under some specific contexts, justifying the term 
entropy functional but it should be stressed that the results in this section are logically independent of the derivation 
of field equations in the previous section. In particular, we need to consider non-null vectors to provide a natural 
interpretation of the extremum value of the functional. 

The term 'on-shell' refers to satisfying the relevant equations of motion, which in this case are given by Eq. (|18[) . 
Manipulating the covariant derivatives in Eq.([l]), we can write 

ski = / d D xv=-g [4v d (p ab cd (VeO e) - 4P ab cd (v d VcO e - T ab ee] 

Jv 

= 4 / d^xVhn, (P ab cd e^cC) + I d D x^—g (2P mb cd R m acd - T ab ) . (29) 

JdV Jv 

In writing the first equality, we have used the condition V dP a bcd — 0. As before, in the first term of second equality, 
n a is the vector field normal to the boundary dV and h is the determinant of the intrinsic metric on dV. (In general, 
the boundary is (D — l)-dimensional. We will soon see that the really interesting case occurs, in fact, when part of 
the boundary dV is null and hence intrinsically (D — 2)-dimensional. This case needs to be handled by a limiting 
procedure and in what follows we will elaborate on the procedure we use.) The second term of the second equality in 
Eq. (f!?9")) vanishes in the absence of matter because, when T ab — 0, the equations of motion reduces to Q a ^ k Rbijk = 0, 
thereby allowing us to interpret the first term as the on-shell value of entropy from the gravity sector. Even in the 
presence of matter, the second term can be expressed in terms of matter variables as an integral over the trace of the 
stress tensor (see the discussion around Eq. (|27[) ) and, of course, is not a surface term. We will, therefore, concentrate 
on the surface term arising from the gravitational sector which reduces to 

SU-sheii = 4 / d D - 1 xVhn a {P abcd ^ b (, d ) — ^ -!- / dP- x xyfhn a (^V^ 6 - ^bC) (30) 

JdV JdV 

where we have manipulated a few indices using the symmetries of p abcd . The second expression after the arrow is the 
result for general relativity; we give this explicitly to show the form of the expression in a familiar setting. Note that, 
when £ a is chosen as the normal to a set of surfaces foliating the spacetime, the integrand has the familiar structure of 
Ui^K + a 1 ) where a 1 = £ Vb£* is the acceleration associated with the vector field £ a and K = — Vh£ b is the trace of 
extrinsic curvature in the standard context. This is the standard surface term which arises in the ADM formulation 
and the cognoscenti will immediately see its connection with entropy of horizons in GR. 

At this stage, we have not put any restriction of the boundary dV or on the choice of the vector field £ a . The 
expression in Eq. (|30p is valid for any vector field £ a — not necessarily null. (Our entropy functional is defined for 
any vector field; to obtain the equation of motion we consider only the variation of null vector fields but having done 
that, we can study the on-shell entropy for any vector field.) The only restriction is that the expression in Eq. (|30|) 
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should be evaluated on a solution to the field equations. It is clear that one cannot say much about the value of this 
expression in such a general context, keeping the boundary and £ a totally arbitrary. Further, even in the case of a 
null vector field £ a , the integrand X in Eq. ([30f changes by T — > f 2 (x)T, under a rescaling £ a — > f(x)£ a which keeps 
the null vector as null. Since the value of the integral can be changed even by such a rescaling, it is clear that a choice 
has to be made for the overall scaling of the null vector field before we can evaluate S[£] on-shell 23]. 

The fact that Eq. (|30j) has no clear interpretation in general should not be surprising since we expect to obtain a 
nontrivial value for the entropy only in specific cases in which the solution has a definite thermodynamic interpretation 
and the surface and the vector field is chosen appropriately. Obviously, making this connection will require choosing 
a particular solution to the field equations, a particular domain of integration for the entropy functional, etc. and 
making other specific assumptions. We shall now calculate the extremum value in specific situations and demonstrate 
that it gives the standard result for the gravitational entropy when the latter is well-defined and understood. We 
will also discuss several features of this issue and, in particular, will demonstrate that in the standard cases with 
horizons, the extremal value of the entropy correctly reproduces the known results, not only in GR but even for 
Lanczos- Lovelock theory. 

The most important case corresponds to solutions with a stationary horizon which can be locally approximated 
as Rindler spacetime. Many of the results which had motivated us to develop the current formalism were proved 
in this specific context and hence this will act as a natural testing ground. In this case, the relevant part of the 
boundary will be a null surface and we will choose £ a to be a spacelike vector which can be interpreted as describing 
the displacement of the horizon normal to itself. To define this properly, we will use a limiting procedure and provide 
the physical motivation for the choice of £° (based on certain locally accelerated observers) thereby leading to a 
meaningful interpretation of the on-shell value of £[£]. 

To set the stage for calculations that follow, we will begin by briefly recalling the notion of Rindler observers in flat 
(Minkowski) spacetime. In Minkowski spacetime with inertial coordinates x l M — (tyi, ^m)' ct = 1, 2, D — l, observers 
undergoing constant acceleration along the direction (the Rindler observers) follow hyperbolic trajectories 0, [l6| 
described by (x^) 2 — (^m) 2 = constant. A natural set of coordinates for these observers is given by x^ — (iR,, N, x±), 
A = 2,3, D — X, where the transformation between x^ and xf^ are given by 

N . , , > i N , , , 
£m = — smmraR.) ; x M = — cosh(ftiR) 

K K 

N = n({x\ A f-{t M f) 1/2 ; ^Itanh- 1 ^) 

x*=xj_ , A = 2,3...,D-1, (31) 

with constant k. The metric in the Rindler coordinates becomes 

ds 2 = —N 2 dt 2 l + dN 2 /n 2 + dh\ , (32) 

where dL\ is the (flat) metric in the transverse spatial directions. It is easy to see that the surface described by 
( x m) 2 ~ (^m) 2 = (or N — in the Rindler coordinates) is simply the null light cone in the <m — plane at the 
origin, and that it acts as a horizon for the observers maintaining N = constant =/= 0. 

In a general curved spacetime, one can introduce a notion of local Rindler frames along similar lines. We first go 
to the local inertial frame (LIF, hereafter) around any event V and introduce the LIF coordinates x\^ = {jbyi, Xyi), 
a = 1,2, D — 1. We then use the transformations in Eq. (f3Tj) to define a local Rindler frame (LRF, hereafter). The 
choice of axis is of course arbitrary and one could have chosen any direction in the LIF as the x\^ axis by a simple 
rotation. In particular, a general null surface Tt in the original spacetime passing through V can be locally mapped to 
the null cone in LIF which — in turn — can be locally identified with the N = surface for the local Rindler frame. 
This local patch 7Ylif C H of the original null surface acts as a horizon for these observers. We will make good use 
of this observation below. (The local nature of the construction is more transparent in the Euclidean description. If 
we choose a LIF around any event and then transform to an LRF, then the null surface in the Minkowski coordinates 
gets mapped to the origin of the Euclidean coordinates. Our constructions in a local region around the origin in the 
Euclidean sector captures the physics near the Rindler horizon in the Minkowski frame.) The crucial fact to notice 
is that locally, Ti is the Killing horizon for a suitable class of Rindler observers. To see this, choose some point 
V € H and erect a D-ad (the D-dimensional generalisation of a tetrad) in the LIF at V , endowed with Minkowski 
coordinates. Let 7Ylif C H denote the part of 7i contained in the LIF. Choose the D-ad in such a way that the only 
non- vanishing components of the generator \ a of 7Ylif are \° an d X • I n other words, with respect to this D-ad, Hlif 
is defined by (x^) 2 — (x^) 2 = 0, where x^ are the Minkowski coordinates in the LIF. Now transform to the local 
Rindler frame using the transformation in Eq. (|31[) and consider the vector v a = (1,0) in the Rindler frame. Clearly 
v a is the Killing vector associated with time translations in the Rindler frame, with norm v a v a = —N 2 , and hence 
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Wlif (given by N = 0) is a Killing horizon for the Rindler observers, generated by v a . (It can be shown that the 
original generator x a of Hlif when transformed to the Rindler frame, is proportional to the Rindler Killing vector v a 
on the horizon Hlif-) 

We will now give a prescription for the evaluation of i5| on _ s heii in a specific LIF (i.e. on 7Ylif C Ti), which extends 
to the entire surface H. in an obvious way. For notational convenience therefore, we will drop the subscript 'LIF' on 
H . Instead of the surface H, consider the surfaces in the local Rindler frame at V given by N = e = constant. Take 
£a _ n a as ^ e un j-j. gpacgiike normal to these surfaces, so that 

n a = C = (0,k,0,0,...); Vh = ey^, (33) 

where a is the metric determinant on the £r = constant, N — constant surfaces. We will evaluate the surface integral 
for S'lon-shoii on a surface with N = e = constant, and take the limit e — > at the end of the calculation. The vector 
displacement field £ a then has the natural interpretation of moving the surface N = e = constant normal to itself and 
previous work has shown that these displacements play a crucial role in the thermodynamic interpretation 0, H[ . In 
our limiting procedure, we use the normal vector to the surface N = e = constant, fix its norm when the surface is 
not null (i.e., e ^ 0) and take the e — > limit right at the end so that it remains properly normalised. In this process, 
we are considering the null surface as a limit of a sequence of timelike surfaces. This is clearly only an ansatz and - 
as we have said before - the final result will be different for a different ansatz. In this sense, it is the end result that 
we obtain which provides further justification for this choice. This prescription can be understood in a natural way 
in the Euclidean continuation (t^ — > Ur) of the Rindler frame, where the N = e surfaces are circles of radius e/ k in 
the t — x plane of the Euclideanised Minkowski coordinates. 

Computing the entropy functional using this vector field, and taking the e — ► limit at the very end, we find that 

S\ n = 47rmc ™ / dD ~ 2x ^ C \^-r> > ( 34 ) 

m=l ^ n 

where x± denotes the transverse coordinates on 7i, a is the determinant of the intrinsic metric on H and we have 
restored a summation over m thereby giving the result for the most general Lanczos-Lovelock case. The proof of 
Eq. (|34[) can be found in the Appendix (| A 2[) . The expression in Eq. (|34[) is precisely the entropy o f a general Killing 
horizon in Lanczos-Lovelock gravity based on the general prescription given by Wald and others [24j and computed 
by several authors [2a |. This result justifies the choice of vector field £ a used to compute the entropy functional 
(as well as the nomenclature 'entropy functional' itself). For a wide class of Killing horizons (Hk) it is possible to 
take the Rindler limit of the near-horizon geometry, and write ds 2 — —N 2 dt 2 + dN 2 / n 2 + dL\ near the horizon, 
where dL\ denotes the line element on the transverse surfaces (and in particular on the N — surface which is the 
horizon; around any point V S 7Yk the transverse directions will be locally flat in the Rindler limit). In this case, k 
is the surface gravity of the horizon (being constant over the entire horizon) and we choose £ a by the above limiting 
prescription. Our construction is then valid over the entire surface 7Yk and the resulting on-shcll value of the entropy 
functional is precisely the standard entropy of the horizon. 

To summarise, we get meaningful results in two cases of importance. First, whenever we have solution to the field 
equations which possesses a stationary horizon with Rindler limit, we have a natural choice for £ a through a limiting 
procedure such that the extremum value of the entropy functional on shell matches with the standard result for the 
entropy of horizon. Second, in any spacetime, if we take a local Rindler frame around any event we will obtain an 
entropy for the locally defined Rindler horizon. In the case of GR, this entropy per unit transverse area is just 1/4 as 
expected. This requires working in a local patch and accepting the notions of local Rindler observers and local Rindler 
horizons about which there is still no universal agreement. (Not everyone is comfortable with deSitter universe having 
an observer dependent entropy let alone Rindler horizon, but we believe this is the correct paradigm.) Finally, in the 
situation in which the boundary dV is nowhere null and £ a is an arbitrary vector field - as we pointed out before - it 
is hard to say anything about the value of S. We hope to investigate this case fully in a later work. 

IV. THE RESULTS IN THE LANGUAGE OF FORMS 

The purpose of this Section is to recast our formalism in the language of forms for the sake of those who find such 
things attractive. This will, hopefully, help in further work because of two reasons. First, the expression for Wald 
entropy [24[ can be expressed in the language of forms nicely. Second, the action for Lanczos-Lovelock gravity can be 
expressed in terms of the wedge product of curvature forms. We shall begin with brief pedagogy to set the stage and 
notation and then will derive the key results. 
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Since the tensor P ab cd has the same algebraic structure as curvature tensor, one can express it in terms of a 2-form 
analogous to the curvature 2-form. We define a 2-form P ab related to our tensor P ab cd by: 

P^Eir^^A^. (35) 

Throughout this section we will assume a coordinate 1-form basis u) % = Ax 1 . If v = e a v a is a vector, then the 
vector- valued 1-form dv is given by 

dv = e a (dv a +uj a b v b ) =e a (X7 b v a )u; b ; u,\ = T a bc u c . (36) 

We will work in the case of pure gravity (that is, T ab = 0) since this is more geometrical and since it does not affect 
the value of the final on-shell entropy functional. Our entropy functional in Eq.JT]) has the integrand (Z?-form): 

l=(4P ab cd V c ^V d e)e, (37) 

where e = (1/ D\)e ai ... aD u> ai A ■ ■ • A u> aD is the natural £>-form on the integration domain V. The first point to note 
is that the £)-form in Eq. (f3"7} is the same as the following: 

X = 4(*P Qb A(dO a A(dO b ) ■ (38) 
where (d£) a = (V c £ a )u; c etc. To see this, we use 

* Pa& = {D \ 2 )\ k\ PabCdecdai - aD -^ ai A " ' A ' (39) 



which allows us to expand the £>-form in Eg. ([38)1 as 

1 = J^^ P ab Cd ^ dai ...a D ^a D ^ a )^a D i b )^ A • • • A (40) 

Since this is a D-form in a D-dimensional space, the right hand side of Eg. (140)) should be expressible as fe where / 
is a scalar. This implies 

2 d 1 

JJj~~2~y ^ >ab € cda 1 ---a D - 2 fflan-i £ )(^a D £, ) = ~j^ f e a 1 - a D (41) 

We now use the identity 

( 6,...6„. ,a,...a i(6] ^ ^ ^ = ( _ 1)S(£) _ j)^...^ ^ (42) 

where s is the number of minus signs in the metric and 5..'. is the alternating tensor with our normalisation. This also 
implies 

e ai - a °e ai ... aD = (-l) s D\. (43) 
Contracting both sides of Eg. (|4"l"T) with e °i-" ^ anc j using the symmetries of e ai ... aD leads to 

/ = 2P afc cd (V,r)(V^ b )^ = 4F Qb cd (V c D(V^ h ) , (44) 

which is exactly what is required in Eo. (|37p . Therefore the entropy functional for any vector field can be written, 
somewhat more geometrically, as 

S[£] = [ 4 (*P Qb A (d£) a A (d£) b ) (45) 
Jv 

We will now do two things. The first is to derive the equations of motion and hence the on-shell value of S[£], and the 
second is to relate this value to the Wald entropy. To derive the equations of motion, it is convenient to introduce some 
new notation which makes the derivation more compact. (Essentially, we will suppress all explicit index occurrences 
in, say Eq. (|4"5|) etc.). The relation to Wald entropy, however, is more easily seen with the indices explicitly in place, 
so we will revert to the explicit notation at that stage. 

To get rid of tagging along the indices, we will first introduce the convention that in a tensor valued p-form T, the 
tensor indices will always be thought of as being superscripts. For example if T is a 2-tensor valued p-form, we have 
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T = e a e T ab where T ab is a p-form and we have suppressed the direct product sign for the basis vectors. We come 
across p-forms like P ab which are antisymmetric in a and b. These can be denoted by the "bi-vector valued" p-form 
P = (l/2!)e a A ebP af> . This will allow us to consistently introduce a dot product for the tensor valued forms. Let us 
now introduce the notation "wedge-dot" A as follows in terms of couple of examples: (1) If A is a 2-tensor valued 
p-form and B is a vector valued g-form then 



A A B = e a (e6 • e c ) A ab A B c = e a A a b A B 



(46) 



That is, the wedge in A acts in the usual fashion on the p-forms and the dot acts on the two nearest basis vectors 
it finds. This last condition also implies that: (2) If A is a bivector valued p-form (A ab = — A ba ) and B is a vector 
valued g-form, then 



AAB = (-1) M+1 BAA, 



(47) 



where the (— l) pq is the usual factor on exchange of the wedge product, and the additional factor of (—1) arises due 
to the dot in A shifting from one index of A to the other. 
In this notation, we have 



* P ab A (d£) a A (d£) h = (-l) D - 2 (d£) a A *P ab A (d£) h = {-\)°- 2 d£ A * P A d£ . 
Further, in this notation one can show that for the bivector valued 2-form P, the following relation holds 



* d * P 



1 

2! 



e a A e fc (V C P 



abc \ , ,d 



(48) 



(49) 



and the condition V C P' 



abed 



is equivalent to d * P = 0. This also means d(*P • £) = ( — 1 



,D-2 



PAd^ (where the 



ordinary dot is defined in the obvious way on the nearest basis vectors), and the entropy functional becomes 



(-i) 



D-2 



d£ A * P A d£ 



d£Ad(*P-€) 



Using the identity 

d (d£ A * P • £) = d 2 £ A * P • £ - d£ A d (*P • £) 
and the definition of the (bivector valued) Riemann curvature 2-form via 



d £ = R •£ = -£• R 



R = ie a A e b R afc 



the entropy functional can be rewritten as 



S[S] 



[ 4 [€ • R A * P • £] - [ 4 (d£ A * P • £) , 

Jv JdV 



(50) 



(51) 



(52) 



(53) 



where dV is the (D — l)-dimensional boundary of the volume V. 

We will now specialize to a subset of vector fields which are null and vary the entropy functional with respect to 
them. The first term in Eq. (|53p is manifestly symmetric in £ and will give a contribution with a factor of 2. (This 
can also be verified using the rules of exchanging the A etc., or by explicitly working with the indices in place.) The 
second term of Eq. (|53"]) will not contribute since the variation 5£ vanishes on the boundary. The condition £ • £ = is 
preserved by adding a Lagrange multiplier term A£ • £ to the entropy functional. The resulting equations of motion 
arc therefore 



-4£-RA *P + A£ = 0, 
which is a set of vector equations. Taking a dot product with £ we get 



(54) 



-4^RA *P-^ = 0, 



(55) 
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which must hold for all null vector fields Eq. ([55"|) can clearly also be written (since d * P = 0) as 

4[d(d£A *P)]-€ = Q ; (56) 

which can also be derived directly by the variation of Eq. ([50|) . It is easy to show that Eq. (|55f is equivalent (after 
bringing back all indices) to the vacuum version of Eq. (fT5|) . To see this, consider 



-4£ ■ R A * P ■ £ = -4£ Q R Qt> A *P b % 
1 



(D-2) 



CRabk^P^^nnk.-kn^ 1 A • • • U) kr> . (57) 



On equating the right hand side of Eq. ()57|1 to fe where / is a scalar and using arguments similar to the ones following 
equation Eq. (j4*0|) , we find 

/ = 2P mb ij R m aij ee = 2P b ijk R a ijk U b , (58) 



which is what appears in the first term on the left hand side of Eq. ([18)) . The rest of the derivation of field equations 
follows as before. 

Having derived the equations of motion by extremising the entropy functional with respect to the null vector fields, 
we go back to the entropy functional for an arbitrary vector field and evaluate it on-shell. Using Eqns. (|5l5|) and (|55[) . 
the on-shell value of the entropy functional for an arbitrary vector field is 

SU_ she ii = -4 / (d£A*P-£). (59) 

J dV 

To show that the on-shell value agrees with Wald entropy for our specific choice described earlier, we take dV to be a 
null surface and denote an arbitrary spacelike cross-section of dV by Ti. Reverting to (semi)index notation, we have, 



5| on _ shcll = -4 / (d£)« A *P°»e 

J dV 



b 

dV 

= ~ 4 L V (v ^ q) (jD^.h pab0dChe °^- aD - 2 ) " k A " ai A " " ■ A ujaD ~ 2 ■ m 

To show that this expression is same as the Wald entropy, we will again employ the limiting procedure used in section 
3; namely, we will consider the Rindler limit of the geometry near the null surface dV. We consider the vector field £ 
to be normal to surfaces of N = constant for the Rindler metric Ea.([3"2")) and take the N — > limit at the end. The 
spacelike cross-section 7i of dV corresponds to the transverse directions of the Rindler metric and has coordinates 
labeled with upper-case indices x^, A = 2, 3, • • • ,D — 1. 

For the vector £ = £ a e a = (0, k, 0, ...), we have Vfc£ a = — (l//t)r^ of which only the component Tq = k 2 N ^ 0. 
The integrand I of Eq. ([6H]) becomes 

I = AN ( 0D~h)! I) p0NCd ^ ai -a D ^° A u> ai A ■ • • A . (61) 

Since the integrand is a (D — l)-form restricted to a surface of constant N, and since the basis 1-form u>° appears 
explicitly, it follows that the remaining basis 1-forms must be intrinsic to the spacelike cross-section Ti.. The indices 
a\ ■ ■ ■ cld-2 hi Eq. (|6ip can then be replaced by A\ ■ ■ ■ Ad-i, and the integrand reduces to 

I = 4N {g m g NN P QN 0N ) V^^-^^V^e Al ... AD _ 2 u: A uj Ai A • • • A u A -^ , (62) 

where e^ 1 ... J 4 D _ 2 ^ s the alternating symbol which takes values +1(— 1) when the indices are an even (odd) permutation 
of 2, 3, • • • , D — 1, and is zero otherwise. The metric coefficients refer to the Rindler metric in Eq. (|32"|) . and cancel the 
N dependence of the integrand. Recognizing the natural (D — 2)-form on Ti. given by 

' M... ifl _ 2 ^A--Aw A ^, (63) 



(D-2)) 
the integrand becomes 



-^"VAe". (64) 
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In a standard derivation of the expression for the Wald entropy, e.g.- Ref. [26|, one would be dealing with the binormal 
to the cross-section TL, defined in terms of two null vectors normal to TL. Since we are dealing with a limiting procedure 
in which the horizon (null surface dV) is approached in the limit N — > 0, we define a 2- form n = (l/2!)n.y u; 4 A u)' J 
which reduces to the standard binormal to TL on the horizon. To do this we use the fact that in the standard case, 
the natural (D — 2)-form e on TL is simply the dual of the 2-form binormal n. We therefore define our 'binormal' via 
the relation 

""- ^e Al ... A „ ^ A '- A " (65) 



(D-2) 



It is now straightforward to check, using \/—g = \ / —gQogNN^/&, that P abcd n a bncd = — ^Pqn° N ■ ^ n Appendix A. 2 



ON 

we have explicitly shown that this quantity is independent of N in the Lanczos-Lovelock case, and hence the TV — > 
limit is trivial. The on-shell value of the entropy functional therefore becomes 

SU-shell = « / P abcd n ab n cd u° A e . (66) 
JdV 

Since a; = dirt, when the on-shell entropy is evaluated in a solution which is stationary, the above integral splits 
into a time integration and an integral over the arbitrary cross-section TL of dV. Then restricting the time integration 
range to from to 2tt/k, we get 

SU-shell = 2^ I P abcd n ab n cd e , (67) 

which is precisely the expression for Wald's entropy when we recall that P abcd = (dC/dR a bcd) for the Lanczos-Lovelock 
type theories (up to a sign which depends on the sign convention used for e) . Using the definition of n»j and following 
the algebra in Appendix A. 2, one can also easily recover the expression in Eq. (|34p thereby proving the equivalence of 
the two approaches. 



V. DISCUSSION 



Since we have described the ideas fairly extensively in the earlier sections as well as the introduction, we will be 
brief in this section and concentrate on the broader picture. 

We take the point of view that the gravitational interaction described through the metric of a smooth spacetime is 
an emergent, long-wavelength phenomenon. Einstein's equations provide the lowest order description of the dynamics 
and one would expect higher order corrections to these equations as we probe the smaller scales. It was shown in 
Ref. 0] that the Einstein equation G® = 87tTq for spherically symmetric spacetimes with horizons can be rewritten in 
terms of thermodynamic variables and is in fact identical to the first law of thermodynamics TdS — dE+PdV. In Ref. 
Q this result was extended to the £q = (1/2)Tq equation for spherically symmetric spacetimes in Lanczos-Lovelock 
gravity where E% was defined in Eq. ([22|) . In fact, the Lorentz invariance of the theory, taken together with the 
equation E® = (1/2)T§ then leads to the full set of equations Ef = (1/2)T£ governing the dynamics of the metric g a b- 
The invariance under Lorentz boosts in a local inertial frame maps to translation along the Rindler time coordinate 
in the local Rindler frame. Hence the validity of local thermodynamic description for all Rindler observers allows 
one to obtain the full set of equations from the time-time component of the equations. Given the key role played by 
horizons in all these, it seems natural that we should have an alternative formulation of the theory in terms of the 
entropy associated with the horizons. This is precisely what has been attempted in this paper. 

The first key result of this paper is an alternative variational principle to obtain not only the Einstein's theory 
but also the more general Lanczos-Lovelock theory. We have shown that there is a natural procedure for obtaining 
the dynamics of the metric (which is now interpreted as a macroscopic variable like the density of a solid) using the 
functional S[£] given in Eq.(TI|. This functional can be defined for any vector field £ a . When we restrict attention to 
null vector fields, and demand that the entropy associated with all the null vectors should be an extremum, we obtain 
a condition on the background geometry that is equivalent to the dynamical equations of the theory. This provides 
an alternative route without the usual problems which arise in the handling of surface terms etc. when the metric is 
varied. Interestingly, our approach selects out the Lanczos-Lovelock type of theories, which are known to have nice 
properties regarding the integrability of the field equations etc. 

As an aside, we want to make a remark on the usual derivations of the field equations in Lanczos-Lovelock theory. 
To illustrate the point consider the familiar D = 4 case. We know that in 4-D, the second order Lanczos-Lovelock term 
£-2 (the Gauss-Bonnet term) is a total divergence, and the higher terms identically vanish. Usually in the literature, 
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one will ignore the Gauss-Bonnet term, since it is a total divergence, and claim that the equations of motion are 
identical to Einstein's equations. One must note, however, that in a situation wherein the Lagrangian contains a total 
derivative term, the conventional action principle is well defined only when all surviving surface terms are held fixed. 
It is well known that in Einstein's GR this can be achieved by adding the Gibbons-Hawking term to cancel certain 
surface contributions. The case of the Gauss-Bonnet term in 4 dimensions proves to be trickier in terms of defining 
consistent boundary conditions. (See Ref. [21( for an attempt to address this issue, and for further references.) So the 
Lanczos-Lovelock theory, in the conventional formalism based on varying an action principle, faces certain difficulties 
associated with the boundary term. Our formalism, of course, reproduces the standard result and the correct entropy 
and — as a bonus — the above difficulties do not arise in our approach since we do not vary the metric to get the 
equations of motion. 

The second result in the paper is related to providing an interpretation of the extremum value of this functional. 
Here we find that one needs to make some ansatz which — though physically well-motivated and natural — is 
logically independent of the derivation of field equations. The ideas we have used are closely connected with previous 
attempts to interpret the radial displacements of horizons as the key to obtaining a thermodynamic interpretation of 
gravitational theories 0, Q ■ When we evaluate the on-shell value of the entropy functional associated with a normal 
to a sequence of timelike surfaces and take the limit when these timelike surfaces approach the horizon, we obtain the 
standard entropy of the horizon. In the case of GR, the result is intuitively obvious since the integrand of the entropy 
functional has the structure nj (K£ J +o J ) which is essentially the surface term in Einstein-Hilbert action from which 
we know that the correct entropy can be obtained. Remarkably enough the same prescription works even in the case 
Lanczos-Lovelock gravity for which there is no simple intuitive interpretation. 

Thus, at the least, we have provided an alternative variational principle to obtain not only Einstein gravity but 
also its closely related extensions, without varying the metric in the functional. This, by itself, is worth further 
study from three points of view. First, it is important to understand why it works. In conventional approaches one 
interprets the extremum value of action in terms of the path integral prescription in which alternative histories are 
explored by a quantum system; here this should correspond to fluctuations of the light cone structure in some sense. 
It is not clear how to make this notion more precise and useful. Second, the matter sector is — as usual — quite 
ugly and nongeometrical and one could even claim that it was added by hand. It is not clear whether the entropy 
functional, including the matter term has a geometric interpretation [191 ]. Finally, the work clearly endows a special 
status to the Lanczos-Lovelock theory as a natural extension of GR within the thermodynamic paradigm. Several 
previous results, especially Ref. [171 ]. have already pointed in this direction. Given the rich geometrical structure of 
the Lanczos-Lovelock theory (compared to, for example, theories based on f(R) Lagrangians), it is worth investigating 
this issue further. 

These results certainly indicate a deep connection between gravity and thermodynamics which goes well beyond 
Einstein's theory. The general Lanczos-Lovelock theory, which is expected to partially account for an effective action 
for gravity in the scmiclassical regime, satisfies the same relations between the dynamics of horizons and thermody- 
namics, as Einstein's GR. This sug gest s that these results have possible consequences concerning a quantum theory of 
gravity as well. In a previous work [13] it was shown that this class of theories exhibits a type of classical 'holography' 
which assumes special significance in the backdrop of current results. 
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APPENDIX A 

In this appendix we provide proofs for some of the results quoted in the text. 

1. Proof of Eq. (Tl9)) 

It is sufficient to prove that if a symmetric tensor S a b satisfies 

s ab ez b = o, (ai) 

for an arbitrary null vector field £°, then it must satisfy S a b — F{x)g a b where F is some scalar function. Consider a 
point V in the spacetime and construct the local inertial frame (LIF) at V , endowed with Minkowski coordinates. The 
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vector field £ a appearing in Eq. (|Aip is arbitrary, and the choice of the D-ad erected in the LIF is also arbitrary up to a 
local Lorentz transformation (LLT) (which has D(D — l)/2 degrees of freedom). We will utilize the (D — 1)(D — 2)/2 
(spatial) rotational degrees of freedom available in the LLT to choose the D-ad axes such that the purely spatial 
components of S ab with respect to these axes, vanish; namely S a p — 0, a, (3 — 1, 2, D — 1. Similarly, we utilize the 
(D — 1) boost degrees of freedom in the LLT to choose the D-ad axes such that So a = 0, a = 1, 2, D — 1. Further, 
since the null vector held £ is arbitrary, we will in turn consider the vector fields given by 

d? a) =S% + 5 a a ,a = l,2,...,D-l, (A2) 

where the components are understood to be with respect to the local Minkowski coordinates. Substituting for this 
choice of £ a in Eq. (|A1|) , we find 

<5oo + 25o Q + S aa = . (A3) 

The middle term drops out because of our choice of D-ad axes. Repeated application of Eq. (|A3[) for all allowed values 
of a then gives us S aa — —Sqq for all a, and combined with So a = = S a p, we obtain 

S ab (V) = F{T) Vab , (A4) 

where F is a scalar depending on the choice of V . Since Eq. (|A4[) is a tensor equation in the LIF, it immediately 
generalises as required, to 

S ab [x) = F(x)g ab (x) , (A5) 



2. Proof of Eq. pi) 

We now consider the algebraic details in the evaluation of the on-shell value of entropy functional. As explained in 
the text, we expect to obtain meaningful result only for certain solutions, when a specific choice is made for the for 
the boundary and the vector field. To do this we define it through a limiting process involving a sequence of timelike 
surfaces and their normals with the limit taken at the end of the calculation. According to the prescription laid down 
in the text, we take £ a = n a to be the unit spacelike normal to the N — e surfaces in the Rindler frame, and taking 
the e — > limit at the very end of the calculation. (This limiting procedure is physically well motivated; in the case of 
standard GR and a Schwarzschild black hole, for example, it will correspond to approaching the r = 2M surface as the 
e — > limit of r = 2M + e sequence of surfaces.) In the Rindler frame, with the metric ds 2 = — N 2 dt^ + dN 2 / k 2 + dL^, 

n o =6, = (0,l/«,0,0,...) ; n a = f - (0, k, 0, 0, ...) ; V^-e^, (A6) 

where h is the metric determinant for the N = e surfaces and a the metric determinant for the = constant, 
£b, = constant surfaces. (When e ^ 0, the surface and its normal are not null but e — is the null Rindler horizon). 
The entropy functional is 

S[C\ = f d^xVhna (4P abc %V b ^ d ) . (A7) 

JdV,e 

Here, d D - 1 x = dt R d D ~ 2 x ± , and V b £ d = -r a bd £ a = -{l/n)T N bd , of which only T N 00 = n 2 e ^ 0. The e on the 
integration symbol reminds us that we are not actually on the given boundary dV, but will approach it as e — > 0. As 
mentioned in the text, our choice of the tensor p abcd is a single tensor ( m ) p abcd — rn^ m 'Q abcd , and appears linearly 
in the expression for S^on-shcii- It is sufficient to analyse this expression using the single ( m )p abcd and to take the 
appropriate linear combination for the general Lanczos-Lovelock case at the end. This will not interfere with the 
process of taking the limit e —> 0. As in the text, we will drop the superscript (m) for notational convenience. The 
integrand for a single m can be evaluated as follows 

— f ^ ( A pNONO ^ff \ 

_ e2 V° f^pN0N0\ 
K ' 

= ^ (-4m 5 W) 

K 

- K\fo (AmQ m m ) . (A8) 
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Consider the quantity Q N0 N0 which — for the ra-th order Lanczos-Lovelock action — is given by 

NO _ ^_ cN0a 3 ...a 2m ( n& 3 i>4 r>b 2m -ib 2m \ I aq\ 

lQ n 2 m °N0b 3 ...b 2m \ It a 3 a i --- ri a 2m - 1 a 2m ) ■ 



' NO 



The presence of both and TV in each row of the alternating tensor forces all other indices to take the values 
2, 3, £)-l. In fact, we have S^-^ = SbIbX'.'A^Z with A i> B i = 2 , 3 > ■■■■> D ~ 1 ( thc remaining combinations of 
Kronecker deltas on expanding out the alternating tensor are all zero since S° A = = 6% and so on). Hence Q N0 N0 
reduces to 

Q N0 N0 = l(^^)si^-i^(R^\..Rf m - 1 f m ) . (AlO) 

In the e — > limit therefore, recalling that R^d \h= ^ D 2 ^Rcd I«j we nn d that Qno N ° ^ s essentially the Lanczos- 
Lovelock Lagrangian of order (m — 1): 



Qno N ° = ^(m-i) \n > (AH) 



and the entropy functional becomes 



S\ H = 2m K [ dt R d D - 2 x ± V^C { ^-_l] . (A12) 

J Ti. 



Restricting the tR integral to the range (0, 2it/k) as usual and using stationarity (which cancels the n dependence of 
the result), we get 

S^\ n - Anm f d D - 2 x ± V^£\°Zl] , (A13) 

where we have restored the superscript (m) in the last expression. Finally, taking the appropriate linear combination 
we find 

K K 

^-2) (A14) 



S\n = ^ Sim) \n - ^mc m J d°- 2 x^C 



(m-1) ' 

This is precisely the entropy in the Lanczos-Lovelock theory. We have also verified that this expression is explicitly 
recovered when working with a spherically symmetric metric of thc form ds 2 = —f(r)dt 2 + dr 2 /f(r) + r 2 dQ 2 D _ 2 - ) 

which admits a horizon at some value of r = ru- The calculation in this case can be done using this metric (instead 
of thc Rindler metric) and considering the limit r — > ru. 

Note that, instead of the above limiting procedure, if we had just foliated the spacetime with null surfaces, chosen 
£ a to be the null normal vector field to thc foliating surfaces and taken the boundary to be one of thc foliating surfaces 
(so that the normal vector n a coincides with the null vector on the boundary) then the surface term S'| on _ s hcii will give 
zero on H. This is transparent in Einstein gravity, where the integrand of becomes <~ n a (— £ a (Vi,£ fc ) + £ h Vh£ a ), 
which vanishes on H where n a and £ a coincide and are null. It can be shown that the same result holds more generally 
and is only dependent on the algebraic symmetries of p abcd . Similarly, if we choose £ a = v a = (1, 0) in the local 
Rindler frame — which is the Rindler time translation Killing vector that becomes null on the horizon — and use 
the same limiting procedure, we get a vanishing entropy. On the other hand, the normalised vector in the timelike 
direction N~ x v a gives the same (correct) result as as our choice. Clearly the final result depends on our choice and 
no general statement can be made. 

Finally, to make contact with results in a more familiar setting, we point out that some of these features are not 
unique to the above context. In fact, a similar caveat also applies to the well known Gibbons-Hawking term which is 
a similar surface term arising in the Einstein- Hilbcrt Lagrangian. Apart from some constant proportionality factors 
which are irrelevant to this discussion, this term is actually arises as the integral of a total derivative (in 4 dimensions) 
as 

A GH - f d A x^jV a (n a X7 b n b ) = - [ d 4 x^V a (n a K) - f d 3 xVh(v a n a )K (A15) 

JV JV JdV 

where n a is the normal to thc foliating surfaces, v a is the normal to the boundary dV of the 4-dimensional region V 
and K = — V^n 6 is the trace of the extrinsic curvature. We now consider the case in which the boundary surfaces are 
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chosen to be members of the set of foliating surfaces (e.g, we can foliate the spacetime by t = constant surfaces and 
choose part of dV to be given by t — t\ and t = t^) so that v a = n a . Then we have 

A GH - f d 3 xVh(n a n a )K ~ f d 3 x^/hK , (A16) 

JdV JdV 

provided dV is not null, so that the normal vector can be assumed to have unit norm n a n a = ±1. This is the familiar 
expression often quoted in the literature. Clearly, the above naive argument breaks down (but the result still holds) 
when the spacetime is foliated by a series of null surfaces (n a n a = 0) and the boundary is one of these surfaces. But 
this case also can be handled by a limiting procedure similar to the one we used for computing our surface integral. 
In fact, our prescription essentially foliates the Rindler limit of the horizon by a series of time like surfaces (like the 
r = 2M + e ^constant surfaces in the Schwarzschild) approaching the null horizon in a particular limit (like e — > 0) . 
In the case of GR, this is equivalent to the standard calculation of integrating the extrinsic curvature (defined by 
this foliation) over the surface and - of course - we get the standard result of entropy density being a quarter of 
transverse area. What is more interesting and nontrivial is that the same prescription works in a much wider context 
and reproduces the Lanczos-Lovelock entropy. 
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